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ABSTRACT 

We present CCD photometry for galaxies around 204 bright [mz < 15.5) Zwicky 
galaxies in the equatorial extension of the APM Galaxy Survey, sampling and area 
over 400 square degrees, which extends 6 hours in right ascension. We fit a best linear 
relation between the Zwicky magnitude system, mz, and the CCD photometry, Bccd, 
by doing a likehood analysis that corrects for Malmquist bias. This fit yields a mean 
scale error in Zwicky of 0.38 mag per magnitude: ie Atti^ ~ (0.62 ± 0.05) ABccd 
and a mean zero point of < Bccd — "Tnz >= —0.35 ± 0.15 mag. The scatter around 
this fit is about 0.4 mag. Correcting the Zwicky magnitude system with the best 
fit model results in a 60% lower normalization and 0.35 mag brighter M* in the 
luminosity function. This brings the CfA2 luminosity function closer to the other low 
rcdshift estimations (eg Stromlo-APM or LCRS). We find a significant positive angular 
correlation of magnitudes and position in the sky at scales smaller than about 5 armin, 
which corresponds to a mean separation of 120ft-~^ Kpc. We also present colours, sizes 
and ellipticities for galaxies in our fields which provides a good local reference for the 
studies of galaxy evolution. 
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1 INTRODUCTION 

Some important aspects of galaxy evolution can only be 
understood by studying the statistical properties of nearby 
field galaxies, in particular its luminosity function (LF). As 
well as providing vital information for galaxy evolution stud- 
ies, an accurate knowledge of the present day LF is needed 
to normalize the number counts of galaxies at fainter mag- 
nitudes, and to understand the clustering and large scale 
structure of the galaxy distribution. Moreover, the colour 
distribution of the nearby galaxies provides a basic refer- 
ence to determine star formation rates in galaxies. 

Much recent work has been directed towards studying 
the evolution of the LF using samples of faint galaxies at 
high-redshift, but a consideration of the ensemble of avail- 
able estimates of the LF at low redshifts suggests that there 
are still large inconsistencies which must be reconciled before 
reliable conclusions can be drawn about the implications of 
deep surveys. 

It is common practice to fit the LF to the so called 
Schechter (1976) form: 



<f>{L) = <!> 



* (i^) ^^p("i^) 



(1) 



were luminosity is related with magnitude in the usual 
way, -^ = ]^q0.4(m,-a/) Determinations of the low red- 
shift B-band LF are available from the Stromlo-APM Sur- 
vey (SAPM, Loveday et al., 1992), the Southern Sky Red- 
shift Survey (Da Costa et al., 1994), the CfA2 redshift sur- 
vey (Marzke et al., 1994), and more recently from the ESO 
Slice Project (Zucca et al., 1997) and the 2dF Galaxy Red- 
shift Survey (Folkes et al., 1999). Determinations in the R- 
band are available from the Las Campanas Redshift Survey 
(LCRS, Lin et al., 1994) and from the Century Survey (CS, 
Geller et al., 1997). The best determinations of the Schechter 
function parameters from each of these surveys are listed in 
Table |I|. The SSRS2 is based on the ESO-Uppsala photom- 
etry of Lauberts & Valentijn (1989), but transformed to the 
B{0) system of Huchra (1976). Da Costa et al. (1994) as- 
cribe the diflerence in the LPs of the two surveys to the 
large colour-term used by Lauberts & Valentijn to relate 
the two systems, and so we therefore use mz for Zwicky 
magnitudes and -B(O) for SSRS2 magnitudes. Da Costa et 
al. (1994) quote the rms diflerence between -B(O) and 6j as 
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Survey 


Band 


Ngal 


z 


M, 


a 


0. 


h^ Mpc-3 


CfA2 


mz 


9063 


0.025 


-18.8 ±0.3 


-1.0 ±0.2 


0.040 ± 0.01 


SSRS2 


B{0) 


2919 


0.025 


-19.5 ±0.08 


-1.2 ±0.07 


0.015 ± 0.003 


SAPM 


bj 


1658 


0.050 


-19.5 ±0.13 


-1.0 ±0.15 


0.014 ±0.002 


ESP 


bj 


3342 


0.1 


-19.6 ±0.07 


-1.2 ±0.06 


0.020 ± 0.004 


2dF 


b} 


5869 


0.14 


-19.7 ±0.06 


-1.3 ± 0.05 


0.017 ±0.002 


LCRS 


r 


18678 


0.1 


-20.3 ±0.02 


-0.7 ±0.05 


0.019 ±0.001 


cs 


Rkc 


1762 


0.06 


-20.7 ±0.2 


-1.2 ±0.2 


0.025 ± 0.006 



Table 1. Luminosity Function parameters derived from currently available surveys. 




Figure 1. Luminosity function estimations in the SAPM and the 
CfA2 as a function of absolute magnitudes M in each magnitude 
system {bj or mz)- The continuous lines enclose the 2-sigma re- 
gion in the SAPM estimation whereas the dashed lines enclose 
the 2-sigma region in the CfA2 estimate. 



< 0.2 mag. We quote the results for the red LFs of the Las 
Campanas and Century surveys here for completeness, and 
do not concern ourselves with the details of the transfor- 
mations between red and blue passbands. We note however 
that Geller et al., 1997 find the LF obtained from the CS to 
be in excellent agreement with a prediction based on colour 
transfmorations applied to the SSRS2 LF. 

Allowing for this adjustment between blue and red mag- 
nitudes, the final six rows of Table hi are all in reasonable 
agreement, with only the CfA2 result showing a large devi- 
ation, particularly in the value of 0*. The difference in the 
LFs deduced from the CfA2 and SAPM surveys is illustrated 
in Figure hi 

At face value, the CfA2 LF seems to have less intrin- 
sically bright galaxies (at least 10 times less M=-21 galax- 
ies). Of course, this depends on the transformation between 
Zwicky mz, APM hj or LCRS R magnitudes. Efstathiou et 
al. (1988) adopted the relation mz = B + 0.29 for the CfAl 
survey (Huchra et al. 1983). This transformation would 
move the CfA2 LF in the correct direction to reconcile it 
with other measurements, although Lin et al. (1994) point 
out that a shift of 0.7 mag in M is required to reconcile 
M* with the other surveys. However, a simple shift in the 



magnitude zero-point would not correct for the apparent dis- 
crepancy in the normalisation, which would suggest that a 
more subtle effect is present in the CfA2 data. 

In this paper we present the results of a photometric 
survey of bright galaxies in the overlap region of CfA2 and 
the equatorial extension of the APM Galaxy Survey (Mad- 
dox et al. 1990a, b; Maddox et al. 1991), which we will use to 
investigate this apparent discrepancy. In Section 2 we begin 
by comparing the APM photometry with the Zwicky mea- 
surements. We describe our CCD observations in Section 3, 
and compare our observations with Zwicky's photometry in 
Section 4. In Section 5 we give some other properties for 
the galaxies in our fields and dicuss the implications for the 
LF. In Section 6 we compare our results with the findings of 
other authors and discuss the possible implications for other 
results. No direct CCD calibration has yet been published 
for this part of the APM Galaxy Survey. We will present a 
detailed comparison of our CCD photometry with the APM 
data down to the survey completeness limit in a separate 
paper. 



A COMPARISON OF ZWICKY AND APM 
MAGNITUDES 



From the information used to calibrate the APM survey we 
would expect B ~ 6^^^ + 0.2, given the mean {B-V) ^ 0.7 
and the colour equation bj — B — 0.28(5 — V) (Blair & 
Gilmore 1982; Maddox et al. 1990b; ahhough this shift could 
increase by as much as 0.07 mag according to the findings 
of Metcalfe, Fong & Shanks, 1995). Adopting the relation 
between mz and B used by Efstathiou, Ellis & Peterson 
(1988), this transformation would imply 



mz ^ bj +0.5, 



(2) 



which is close to the shift in M, adopted by Lin et al. (1996), 
and, if taken as they stand, would suggest that the results 
of the two surveys might be consistent, at least in A/*. 

We investigated the usefulness of this relation by con- 
sidering the subset of ~ 100 galaxies found in the part of 
the APM Galaxy Survey equatorial extension which overlaps 
with the Southernmost region of CfA2 (The S+3 sample of 
Marzke et al., 1994). This sample is drawn entirely from vol- 
ume V of the Zwicky catalogue. We have used the publically 
available CfAl catalogue as our source for Zwicky galaxies, 
with magnitudes corrected as described in CfAl. The CfAl 
includes redshifts only for galaxies with mz < 14.5, but has 
magnitudes and positions for mz < 15.5. Marzke et al find 
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15.5 



Figure 2. Comparison of ttie apparent magnitudes for a sample 
of ~ 100 galaxies in the overlap region of the SAPM and the 
CfA2. Magnitudes in the SAPM catalogue correspond to APM bj 
estimations while magnitudes in the CfA2 are Zwicky magnitudes, 
mz- The continuous line corresponds to bj = m^ and the dashed 
line is the best linear fit, mz = bj — 0.6 it 0.5, showing both 
an offset and a large scatter. We would expect to find at least 
mz = bj + 0.5 



that an improvement in the quality of the plate material 
used for the more recent plates of the equatorial extension 
coupled with changes in the plate copying process could re- 
sult in a small change in the saturation correction required 
for bright galaxies, and that this could produce an effect in 
this direction (Maddox, private communication). However, 
account was taken of such effects in the construction of the 
equatorial extension, and any residual effect is likely to be 
much smaller than the discrepancy shown here. 

The transformation of Zwicky magnitudes to B deduced 
by Efstathiou, Ellis & Peterson (1988) are based on a com- 
parison of the CfAl survey data with the Durham-AAT 
redshift survey using 139 galaxies (Peterson et al. 1986). 
These authors find that the different volumes of the Zwicky 
catalogue have large variations in the zero-points, although 
volume V is found to be representative of the calibration 
as estimated from all volumes. These data are limited to 
mz < 14.5, and the luminosity function parameters inferred 
for the CfAl are consistent with the parameters listed for the 
other surveys in Table 1. It is not possible to draw a direct 
comparison of the APM data with Zwicky data brighter than 
mz ~ 14, due to heavy saturation of galaxies this bright in 
the APM data. The most plausible inference from the com- 
parison described above would therefore be a difference in 
the magnitude scales of the two surveys for 6j > 14, with 
the calibration of Efstathiou, Ellis & Peterson (1988) being 
appropriate for the Zwicky system at brighter magnitudes. 



this sample to be representative of the whole of the South- 
ern part of CfA2 (their Figure 3). We inspected this sample 
of galaxies on the film copies of the UKST Illa-J Sky Sur- 
vey plates, and examined the image maps as reconstructed 
from the APM survey data. We removed from our sample 
all those galaxies which had either been broken up into mul- 
tiple image components or which were composed of multiple 
components which had been merged into a single object by 
the APM software. The distribution of these objects in the 
m,z,b] plane is shown in Figure eI A simple least squares fit 
to these data with the slope constrained to be unity gives a 
zero-point shift of 



mz = b3 -0.6 ±0.5, 



(3) 



which is more than a magnitude in the opposite sense to 
that implied by equation H As will be noted below a detailed 
calibration requires a correction for Malmquist bias. 

The APM Survey is internally calibrated by matching 
images in plate overlaps, with the overall zero-point fixed by 
a number of CCD frames distributed over the survey region 
(Maddox et al., 1990b). At bright magnitudes (bj < 17) the 
calibrations are less well determined due to a combination 
of variations in galaxy surface brightness profiles and the 
smaller number of bright galaxies found in the plate overlap 
regions. The equatorial survey extension has been calibrated 
by matching to the original survey using plate overlaps and 
the adopted zero-point for the whole survey taken from the 
original CCD photometry. Maddox et al. (1990b) also used 
their CCD frames to perform an internal consistency check 
on the quality of the plate-plate matching technique, and 
found the residual zero point errors on individual plates to 
be < 0.04 mag. The individual uncertainties in galaxy mag- 
nitudes in the APM system are therefore expected to be 
much smaller than effect shown in Figure Q. It is possible 



3 CCD DATA 



3.1 Observations 



We used the galaxies from the sample discussed in the previ- 
ous section as the basis for a CCD survey to investigate the 
above discrepancy by providing an independent calibration 
for both surveys. This overlap region is essentially defined as 
211^50 < Q < 3'MO and -0.25° <5 < 0.25°. We obtained im- 
ages with the 2.5m Isaac Newton Telescope(INT) and 1.0m 
Jacobus Kapteyn Telescope(JKT) in October 1994. The de- 
cision to use two telescopes was motivated by the number 
of close groups of Zwicky galaxies in the region, so that we 
use the 10' field of view available at the INT to image clus- 
ter fields containing groups of bright galaxies and the 6' 
field of view of the JKT to image individual galaxies. We 
used identical Tektronix 1024 x 1024 detectors (TEK3 and 
TEK4) and Harris B and R filters on the two telescopes. 
We obtained observations of 58 fields in two nights at the 
INT and 73 fields in three nights at the JKT with exposure 
times of 360s and 600s, respectively. We observed a number 
of photometric standards from the list of Landolt (1992), at 
hourly intervals throughout each night, as well as the field 
of the Active Galaxy AKN120 which contains a number of 
photometric standards (Hamuy & Maza, 1989). 

The data were reduced using standard techniques as im- 
plemented in the IRAF ccdred packages, with the exception 
that we used a modified version of the overscan correction 
routine to compensate for a saturation of the preamplifier 
used with TEK4 on the JKT. This effect manifested itself as 
a sudden drop in the overscan level following readout of par- 
ticularly bright stellar objects in the field, and subsequent 
exponential recovery to the normal level as the next ~ 100 



© 0000 RAS, MNRAS 000, 000-000 



4 E. Gaztanaga & G.B. Dalton 






30 



-.. 18 



m 16 



_ 14 
o 

8 
I -0.2 



^ -0.4 



oa 



"T 1 1 1 1 1 1 1 1 1 r 




.^'■.r 



s?- 



J' 



H — \ — \ — — \ — \ — \ — — \ — \ — h 



JlxIpitfltC 



r 



± 



J I I I I I L_ 



± 



14 



16 18 

Bdit (24.0) 



20 



Figure 3. Comparison of total magnitudes BjfiT for different 
isophotos: / = 26.1 mag arcsec"^ to / = 24.0 mag arcsec"^. 



rows were read out. We found that this problem could be 
corrected by fitting the overscan regions on either side of the 
drop for those fields where the effect occurred. 

Extinction coefficients were derived each night, giving 
values in the range 0.10 < ks < 0.12. We determined zero- 
points for the two combinations of telescope and detector to 
be 



Bo, INT = 24.16 ±0.02, 
and 



Bo 



21.77 ±0.02. 



(4) 



(5) 



We also obtained i?-band images of our standards and tar- 
get fields, with zero points determined to be Ro.int ~ 
24.33 ± 0.1 and Ro,jkt = 22.11 ± 0.02, but we were unable 
to determine any significant colour term from our standard 
stars, consistent with previous experience with this combi- 
nation of CCD and filters which are designed to be very close 
to the Johnson-Cousins system. 



3.2 Image Detection and Photometry 

We used the STARLINK PISA image detection and photom- 
etry software for our photometric analysis. This is essen- 
tially the same as the image detection software used in the 
construction of the APM Galaxy Survey (Irwin, 1986). The 
basic input parameters for the detection are the threshold 
intensity per pixel or surface brightness. Is, and the mini- 
mum size of the object to be detected. As- Amongst other 
parameters PISA returns total area At, the ellipticities and 
the isophotal magnitude Bi or the corresponding total mag- 
nitude B resulting from a curve of growth analysis for each 
detected object after removal of any overlapping objects. 



INTl 
INT2 
INT3 
INT4 
INT5 
INT6 



; arcsec 

26.1 
25.7 
25.4 
25.0 
24.6 
24.4 



arcsec 
96 
34 
34 
17 
10 
10 



JKTl 


25.0 


12 


JKT2 


24.5 


12 


JKT3 


24.1 


6 


JKT4 


23.7 


6 


JKT5 


22.9 


3 



Table 2. Input parameter sets used for image analysis. 



Given that our CCD survey was designed to provide 
a calibration of both the APM and Zwicky data, we were 
necessarily interested in a wide range of magnitudes (14 < 
B < 20) and image sizes. For this range of objects there 
was no unique combination of Is and As that could deblend 
the faint objects without breaking the bright ones. In or- 
der to automate this process as much as possible for the 
whole range of magnitudes, we ran PISA several times with 
different combinations of Is/ As. These are listed in Table H. 

We chose a larger (smaller) area for the fainter 
(brighter) isophotes so as to select similar objects in all runs. 
Objects in the final catalogue were selected from the INT4 
and JKT4 runs. The information from different isophotes 
was then used to perform an automatic rejection of broken or 
contaminated images. After rejection, the total magnitude 
and size were the largest remaining estimates of B and At, 
which typically correspond to the faintest isophote left. The 
error in B was taken to be the variance in the different esti- 
mates for total magnitudes. Objects with rejected isophotes 
were automatically labeled and checked visually. We refer to 
total magnitudes estimated in this way as Bimt and Bjkt 
for the INT and JKT sets, respectively, or Bccd in general. 
These effectively correspond to total magnitudes detected 
at isophotes 26.1 mag arcsec"^ and 25.0 mag arcsec"^, re- 
spectively, unless stated otherwise. 

To check if our isophotes are low enough we have 
also computed total magnitudes determined by PISA us- 
ing higher isophotes. We find a small zero-point shift of the 
total Bimt scale as a function of the isophote I, but for a 
given / this shift is not a function of B over the wide mag- 
nitude range considered here. This is illustrated in Figure H 
which shows a change of « 0.2 mag in the mean Bint{I) 
when we change the isophote from 7 = 26.1 mag arcsec"'^ 
to 7 = 24.0 mag arcsec"'^. Changing the isophote from 
7 — 26.1 mag arcsec"'^ to 7 = 25 mag arcsec"^ introduces 
a change of only 0.1 mag. Changing the detection isophote 
from from 7 = 26.1 mag arcsec"^ to 7 = 24.4 mag arcsec"'^ 
also introduces a change in Bimt of about 0.2 mag. We 
therefore conclude that the residuals associated with our fi- 
nal choice of detection isophote are small compared to the 
shifts illustrated in Figure 0. 

This analysis also implies that there should be a mean 
zero-point shift in the JKT data relative to the INT data of 
0.1 mag due to difference in the isophotes used. We therefore 
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Figure 4. The top panel shows the red magnitudes RcCD for a 
subsample of INT galaxies as a function of the blue Bcq£, ones. 
The bottom panel shows the colour BccD-RcCD evolution as a 
function Bqcd- 



apply this shift to all objects observed with the JKT to 
transform these data to our B26.1 system. 

Finally, we investigated the possibility that there could 
be a small effect due to the change in mean redshift of the 
samples at fainter magnitudes by searching for a change in 
the mean galaxy colour. Figure H shows (Bint-Rint) as 
a function of Bint- There is only a small colour evolution 
within the errors. A linear fit to the binned data in Figure H 
yields B - R = 0.016 B + 1.23 with a mean B ~ R= 1.50 
(the mean weighted by each galaxies is B — R = 1.61 as is 
dominated by the faint objects). 



4 COMPARISON WITH THE ZWICKY 
CATALOGUE 

In Figure H we show the number counts histograms for the 
204 Zwicky galaxies in our sample to the different magnitude 
systems: Zwicky (top), APM (middle), and CCD (bottom). 
Comparison of the APM and CCD data shown here suggests 
that the effect shown in Figure His unlikely to be an artefact 
of saturation effects in the equatorial APM Survey data at 
bright magnitudes (see Section 0) . 

We searched for possible correlations between the mag- 
nitude differences, Been — mz, and other properties of our 
sample. The top panel of Figure o shows the Been — rnz 
difference as a function of CCD colour B — R. The colour 
of Zwicky galaxies is similar to the mean colour in Figure 
0, B — i? ~ 1.5 and there is no apparent correlation with 
BccD — mz within the scatter. This is a useful check, as any 
large differences that were due to processing errors might be 
expected to show up as objects of extreme colour. The bot- 
tom panel of Figure pi shows the scatter in Bccd — mz as a 
function of right ascension. Objects were observed in order 
of increasing RA on each night of our observing run, and so 
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Figure 5. Histograms comparing the numbers of galaxies as a 
function of apparent magnitude B. Each panel shows in black, 
from top to bottom, the Zwicky, APM and CCD magnitudes for 
the same objects. In the lower two panels we also show the Zwicky 
histogram as a background for comparison. 
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Figure 6. Bqcd ~ "^z as a function of CCD colour B — R (top 
panel) and right ascention (bottom). 



we would expect temporal drifts to show up as a correlation 
with RA. Again, there is no evidence of any trend. We also 
investigated possible variations with Galactic latitude (not 
shown here), and again found no evidence for trends in our 
data. 

The top panel of Figure shows the B — mz error as 
a function of the mean surface brightness for all Zwicky 
galaxies. Surface brightness is defined here as the ratio 
of isophotal magnitudes to the area above a threshold of 
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Figure 7. Bqcd — "mz as a function of surface brightness (top 
panel) and run number or time of observation (bottom). 



25.0 mag arcsec"'^ in the INT (triangles) and with a thresh- 
old of 24.1 mag arcsec"^ in the JKT (circles). The difference 
in the threshold explains the systematic shift in the bulk 
value of mean high surface brightness of the JKT images 
(as we are closer to the galaxy core). As can be seen in the 
Figure, there is no apparent variation of magnitude error 
with surface brightness. 

The bottom panel of Figure M shows the B — mz er- 
ror as a function of the sequential run number. This is close 
to observational time for the INT (triangles) or JKT (cir- 
cles) when considered separately (in practice, some of the 
INT and JKT observations were made on the same night). 
There is no apparent variation of magnitude error with run 
number. 

Figure ra is the main result of this paper. It shows the 
Zwicky magnitude error Bccd — rnz as a function of Bccd, 
for all Zwicky galaxies with Bccd < 17.5. For the Zwicky 
magnitudes we have used a constant error: Amz = 0.05 mag 
(as Zwicky quoted magnitudes with a precision of 0.1 m ag). 
For the CCD data we used the error described in Section 3.2 



A very similar trend is found for the separate INT (closed 
circles) and JKT images (open squares). 

The dashed line in Figure H shows a direct least square 
fit to the data: 



Bccd — mz ^ 0.6 (Bccd — 15.1) 



(6) 



which has a scatter of 0.55 mag. 

This fit is the result of the interplay between the scatter 
in the two magnitude systems and the Zwicky survey limit 
(shown as dotted line in Figure H). As a result of this scatter, 
objects with faint Bccd and bright mz can be included in 
the survey, but there is a deficit of objects with faint m.z 
and bright Bccd- i.e. the fit suffers from a Malmquist type 
of bias. However, for a linear relation, it is apparent that 
Malmquist bias is insufficient to account for all of the effect 
shown in Figure H. We next develop a scheme to correct this 
fit for the effects of Malmquist bias. 
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The dashed line corresponds to a direct least square fit to the 
data. The continuous line shows the fit corrected for Malmquist 
bias (see text). 



4.1 Malmquist bias correction 

Different magnitude system are subject to different system- 
atic errors, and even in the best of the situations intrinsic 
differences in the morphology, environment and spectrum of 
the galaxies tend to introduce stochastic fluctuations in any 
magnitude system. We want to find a best fit linear relation 
between the Zwicky, mz, and CCD, Bccd, system: 

mz = \BccD + Z, 



(7) 

where A will account for any scale difference and Z is a zero 
point shift. In general, both mz and Bccd are stochastic 
variables and equation H is just a mean relation. As is com- 
mon practice we will assume that there is Gaussian scatter 
around the above mean relation (due to the accumulation 
of multiple uncorrelated factors). That is, given a measured 
magnitude mz, the error is given by: 



P(A) 



iV 



exp 



A^ 

'2cr2 



(8) 



where A = m,z — rnz is a stochastic variate, rnz is some 
mean best fit value in the linear relation of equation M, a 
is the rms error and A*' is a normalization factor. For a 
sample that is not magnitude limited we have A'^ = \/2Tva. 
For a magnitude limited sample, where mz < rn%"^ , the 
probability is the same but has a different normalization: 



N 



dmz exp 



'2cr2 



(9) 



because not all magnitudes are possible, so that A = 
TTiz — mz < 771^"" — rnz- Thus, we can write P(mz) in terms 
of the complementary error function: 
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P(A) = 


2 exp 






-\/27r a erfc 


•J2a 



(10) 



result. 

We are now able to perform a likelihood analysis to 
find the best fit values of A and Z in the linear relation of 
equation M. We define a likelihood as: 



£ = J]P(A,) 



where i runs over all galaxies in the survey, and: 
A, =mlj -{XBhcD + Z) 



(11) 



(12) 



with m''z and Bq^q the measured Zwicky and CCD magni- 
tudes for galaxy i. In analogy with the standard x^ t^st we 
define a Malmquist bias "corrected chi-square": 



2 

XMaln 



E 



^ + 21og(erfc^ 



+ Z- 



m}^" ^ 



^/2g 



(13) 



Note that the measured magnitude errors, Oi, are added in 
quadrature to the stochastic error in the linear relation, o, 
which is one of the parameters we want to determine with 
this fit. The best fit values that we find on maximising the 
likelihood, using the data in Figure |8|, are: 

A = 0.62 ±0.05 (14) 

Z = 5.9 ±0.7 
G = 0.40 ±0.05, 

the errors correspond to approximate 99% confidence levels. 
The values in each uncertainty range are strongly correlated: 
Lower values of A (eg larger scale errors) are (linearly) cor- 
related with larger values of Z, and both are obtained for 
the smaller a. The inverse relation gives: 

BccD ^ l.eimz - 9.5. (15) 

Figure H shows as a continuous line the best fit model for 
the corresponding best fit magnitude difference: 

BccD -mz = (0.38 ± 0.02)(Bccc - 15.53) ± 0.4 (16) 

Figure ^ shows the residuals of this fit. For comparison 
we display a Gaussian with same dispersion, a = 0.4. As 
mentioned above the Gaussian does not provide a good fit 
because of the Malmquist bias, which produces a deficit of 
faint objects. It is not possible to show in this Figure a 
comparison with the Malmquist bias corrected version for 
the error probability of equation nol because this depends 
not only on the differences. A, but also on the measured 
value, B^cD- 

Thus, after correcting for the effects of Malmquist bias, 
the above analysis indicates both a zero-point shift and a 
change of magnitude scale. The zero point different Zo could 
be obtained by taking the mean of the magnitude differences 
over raz magnitudes (which are the ones that define the 
survey limit): 



Zo =< BccD — mz >= 



-0.35 ±0.15. 



(17) 



which is in good agreement with Efstathiou et al. (1988). 
The scaling relation between the two magnitude systems is 
then: 






0.15 



0.1 - 



0.05 




Figure 9. The distribution of magnitude errors from the mean 
linear fit compared to a Gaussian with dispersion a = 0.4. The 
Gaussian does not fit the distribution because of the Malmquist 
bias, which produces a deficit of faint objects. 



A 



Amz 



_ 0.62 ± 0.05. (18) 

ABccD 

We can also test the above model for the scale error 
by estimating the magnitude correlations as a function of 
projected separation. Figure hd shows the mean correla- 
tion < A{9i)A{9j) > as a function of the pair separation 
\9j — 9i\, in arcmin. The uppermost panel shows (as con- 
tinuous lines) the autocorrelation for magnitude differences 
in the Zwicky system, Af — m^ — rnz, where niz is the 
Zwicky magnitude for the ith galaxy and rriz =< mz > 
is the mean Zwicky magnitude in the Survey. The middle 
panel shows the corresponding autocorrelation for the CCD 
system: Ai = Bqcd ~ Bccd- The lower panel shows the 
autocorrelation of the magnitude errors: e^ = Bqqj^ -~ rn^z- 
The long-dashed in each case show the zero values for ref- 
erence. The short-dashed line in the lower panel shows the 
prediction based on applying equation hq to the data, which 
is extremely close to the observed result. 

We can see in Figure hO that there is a significant an- 
gular correlation between nearby magnitudes and positions 
with 9 < 5' . This correlation is followed by the errors, indi- 
cating that the above model is valid. The Zwicky system 
shows smaller correlations and smaller variance than the 
CCD system. This can also be understood in the context 
of the model, as in the Zwicky system the "effective" mag- 
nitude scale is about a factor of two smaller (equation hsl). 

Note that at the typical depth of the survey, 2? ~ 
SOh^^ Mpc, the above magnitude correlations are only sig- 
nificant at physical scales smaller than < lOO/i^^Kpc. 
This correlation has little effect on typical galaxy cluster- 
ing scales, ro — 5/i~^Mpc, but might be relevant for the 
inversion of angular clustering on smaller scales. For ex- 
ample, Szapudi & Gaztafiaga (1998) find that on scales 
6 < 0.1° there is a significant disagreement between the 
APM and the EDSGC that is attributed to differences in 
the construction of the surveys, most likely the dissimilar 
deblending of crowded fields. At the depth of these Surveys, 
V ~ 400 /i~^ Mpc, the above magnitude correlations corre- 
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Figure 10. Mean correlations in magnitude diflFerences < 
A{6i)A{6j) > for pairs of galaxies separated by distance \6j —6i\, 
in arcmin shown for Zwicky magnitudes (top panel), CCD mag- 
nitudes (middle panel), and for the difference e between Zwicky 
and CCD magnitudes (lower panel). The short-dashed line shows 
the prediction using the scale error model e = 0.6AB. The long- 
dashed line shows zero correlation for reference. 



spond to ^ < 1' and could also have a significant efi^ect on 
the angular clustering and its interpretation on these small 
scales. 



4.2 An illustration 

An illustration of the scale error in the Zwicky system can 
be seem in the galaxies shown in Figures 02 and hll. Table 
ra gives the Zwicky and CCD magnitudes for each of the 
Zwicky galaxies as labelled in the Figures. As can be seen 
from the table, the range of Zwicky galaxies in Figure O 
is Amz = 0.3 which is almost 6 times smaller than the 
CCD range ABqcd = 1-7. The range for the whole clus- 
ter is Amz ~ 1.3, almost a factor of 3 smaller that the 
CCD range: ABccd = 3.5. These illustrations are strongly 
suggestive of an observer-bias effect, whereby some fainter 
galaxies are included in the same due to their proximity to 
brighter objects, e.g. object #9 in Figure hi. 



5 THE PROPERTIES OF NEARBY FIELD 
GALAXIES 

We will now present some further properties of the galaxies 
in our sample: colours, sizes and ellipticities, and discuss the 
implications for the Luminosity function. This will help us 
to understand the systematic effects present in the Zwicky 
magnitude system. As mentioned in § 1, these local prop- 
erties are interesting in the context of galaxy evolution and 
star formation rates. 



Figure 11. A CCD image (in B) containing several Zwicky galax- 
ies (labeled with numbers) in a rich cluster. 



Galaxy # Zwicky mz CCD B 



1 


15.0 


15.0 


2 


15.0 


16.1 


3 


14.8 


15.1 


4 


14.7 


14.4 


5 


14.8 


14.9 


6 


14.9 


14.8 


7 


14.0 


13.2 


8 


15.2 


14.9 


9 


15.3 


16.7 



Table 3. Comparison of CCD and Zwicky magnitude for galaxies, 
as labeled in Figures hll and ha. 



There are still only few samples that are both homoge- 
neous and large enough to provide estimates of the statisti- 
cal properties of nearby samples, most of which have been 
selected from photographic plates. Besides the redshift cata- 
logues listed in Table 1, one of the more extensive catalogues 
of bright galaxies is contained in the Second Reference Cat- 
alogue of Bright Galaxies (de Vaucouleurs, de Vaucouleurs 
& Corwin 1976, known as RC2), which gives magnitudes, 
colours, ellipticity and morphology for well over a thousand 
galaxies to a limiting isophote of around 25.0 mag arcsec"^. 
The problem with this catalogue is that it is a mere compi- 
lation of data and was not intended to be complete to any 
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Figure 12. A CCD image (in B) containing several Zwicky galax- 
ies (labeled with numbers) in a rich cluster. 



20 - 



m, < 15.5 




1000 
Area (arcsec^) 

Figure 13. Frequency distribution (in per cent) of galaxies as a 
function of their area (in arcsec^ ) for galaxies selected in differ- 
ent ways: (a) Zwicky sub-sample (top), (b) bright galaxies with 
CCD magnitudes i? < 16 (middle top) and (c) faint galaxies with 
CCD magnitude 19 < S < 20 (middle bottom), (d) area versus 
ellipticitios for all CCD magnitudes 14 < B < 20 (bottom). 



specified limiting magnitude, diameter, or redshift. More- 
over it is based on photographic plates. 

Our sample is homogeneous enough and extends over a 
large enough area ( ~ 400 square degrees) to provide a fair 
sample. Thus the new results based on the CCD magnitudes 
should be a good local reference for the fainter studies of 
galaxy evolution. 

We will compare the properties of galaxies in the Zwicky 
{mz < 15.5) sample with the corresponding properties of 
galaxies in the INT fields selected using the CCD blue mag- 
nitude B = BccD- We choose two magnitude cuts: B < 16 
which includes most Zwicky galaxies and 19 < B < 20, 
which includes faint galaxies in the same fields. 



5.1 Sizes and Ellipticities 

Figure [13 shows an histogram of the frequency distribution 
of galaxies as a function of its area (number of galaxy pixels 
above detection threshold in the CCD image). Top panel 
shows Zwicky selected galaxies only. Middle and bottom 
panels include all the galaxies in the INT fields selected with 
CCD magnitudes of B < 16 and 19 < B < 20, respectively. 

As can be seem in the Figure, Zwicky selected galaxies 
have a long tail of small objects which is not present in the 
bright subsample of B < 16 CCD selected objects. This 
again illustrates the scale error (and selection biases) in the 
Zwicky system mentioned above (see also § 6). 

The distribution of sizes in logarithmic scale for both 
CCD sub-samples are well approximated by Gaussians 
(shown as continuous lines in the figure). The mean size 
is about 60 arcsec^ for 19 < B < 20 and 1600 arcsec'^ for 
B > 16, with a rms dispersion of about 20% and 27% re- 
spectively (recall that these areas correspond to a nominal 



threshold of 26.1 mag arcsec" , fainter thresholds could be 
needed to sample the size of the lower surphace halos) . 

Figure Pl3 shows the corresponding frequency distribu- 
tion for the ellipticities as measure in the galaxy shapes 
(with a threshold of 26.1 mag arcsec"^). These Figures are 
in rough agreement with the results by Binney & Vau- 
couleurs (1982) over the Second Reference Catalogue (RC2). 

The local distribution seems remarkably similar to 
the faint one, given the large differences in sizes shown 
in Fig.h3. This is a good indication that our isophote 
(26.1 mag arcsec^'^) is low enough, as otherwise we would 
expect a large excess of round faint objects. There seems to 
be nevertheless a slight (~ 5%) excess of round faint objects. 
This is probably not due to the seeing (or pixel resolution) 
as most of the faint objects here have more than 30 arcsec'^ 
(or 100 pixels) of area (see Fig.|l3|). The lack of correlation 
between ellipticities e and areas A is illustrated in the bot- 
tom panel of Figure llSl A least-square-fit to the points give 
e~ 0.277 + 0.018 logio(^). 



5.2 Colours 

Both the morphology and a more detailed study of colours 
will be presented elsewhere. Here we just show the colour 
frequencies and discuss some of its implications. 

Figure [13 shows the frequency distribution of CCD 
colour (B — R) for all Zwicky galaxies (top) in our sample. 
We then show the corresponding distributions for the same 
Zwicky galaxies separated into two sets, corresponding to 
the JKT (middle) and INT (bottom) CCD frames. As men- 
tioned above, the INT field of view (10') was larger than the 
JKT (6') and so we used the INT to target groups and clus- 
ters of galaxies (which could take up several overlaping CCD 
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mz < 15.5 (All) 




B-R 



Figure 14. Frequency distribution (in per cent) of galaxies as a 
function of their ellipticity for galaxies selected in different ways: 
(a) Zwicky sub-sample (top), (b) bright galaxies with CCD mag- 
nitudes B < 16 (middle) and (c) faint galaxies with CCD magni- 
tude 19 < B < 20 (bottom). 



frames), while the JKT was used for more isolated galaxies. 
As can be seem in the Figure the INT galaxies are signifi- 
cantly redder, B- R~ 1.57, than the JKT, B - R~ 1.36. 
There are about ~ 100 objects in each subset which roughly 
corresponds to a random sampling of the total 600 Zwicky 
galaxies in our survey area (the 600 targets were split in 
half between the two fields and then more or less randomly 
selected during each night). The mean colour of our 200 
Zwicky galaxies is B — R ~ 1.47. Neither the mean colour 
or the frequency distribution change significantly when we 
exclude all the galaxies that are in the same CCD frame 
(eg to avoid nearby ellipticals). This indicate that selection 
effects are not important for this distribution. 

Given the large area covered (several hundreds of square 
degree) , these results should be a good estimate of the over- 
all mean colour of bright local galaxies. How do these val- 
ues compare with previous studies? This is a difficult ques- 
tion because it requieres both accurate colours and the ac- 
curate fraction of galaxies in different enviroments (eg dif- 
ferent morphological types). Previous studies were limited 
to inhomogeneous samples (eg RC2) or to small numbers of 
galaxies. For example Kennicutt (1992) presented results for 
8 early-type galaxies and 17 spirals to irregulars. Our mean 
B — R ~ 1.47 can be compared againts the synthetic B — R 
colours presented by Fukugita et al. (1995). We have used 
Harris filters which were designed to be very close to the 
Johnson-Cousins (with R Johnson) system. In these broad- 
bands. Table 3 of Fukugita et al. show Br = 1.67 for El- 
lipticals, 1.48 for SO, 1.2 - 1.1 for Spirals and 0.6 for Ir- 
regulars. The range seems in rough agreement with Figure 
na, although we find a significant number of objects with 
B — R > 2. Notice that these objects are mostly in the INT 
sample, that is in clusters or groups. While most of the bluer 
objects with B — R < 1 are in the JKT fields, that is, they are 
more isolated (by angular separations larger than 6', which 
corresponds to a mean separation larger than 140h^^ Kpc. 
Nevertheless, this is seems to be a small effect. 



Figure 15. Frequency distribution (in per cent) of galaxies as 
a function of CCD colour B — R for Zwicky galaxies selected in 
different ways: (a) All Zwicky in our sample (top), (b) Zwicky 
galaxies in our JKT sub-sample (middle) (c) Zwicky galaxies in 
our INT sub-sample (bottom). The same Gaussian distribution 
is plotted as a solid line in each panel 



Note that Fukugita et al synthetic colours seem to be 
slightly less red than the Kennicutt (1992) observations (ac- 
cording to Table 2 in Fukugita et al., they are about 0.1 
redder in B-V and this could be larger in B-R). Also no- 
tice that we are using total magnitudes, while the synthetic 
colours are based in small aperture spectra. Galaxies could 
have quite a different colour distribution (or spectra) in their 
low-surface halos. 

Figure hq shows the frequency distribution of CCD 
colour (B—R) for the Zwicky galaxies (top panel) in compar- 
ison with the same galaxies selected with CCD magnitudes: 
with B < 16 (middle panel) and a sub-set of INT galaxies 
with 19 < B < 20 (bottom panel), which corresponds to the 
points in Figure m The continuous line in all cases shows for 
reference a Gaussian distribution with mean B — R = 1.45 
and rms deviation of a = 0.4, which roughly matches the 
Zwicky frequencies. 

It is interesting to notice the peak of local galaxies with 
B — R = 1.5 and the spread and shift to the red of the faint 
objects in the bottom panel. These magnitudes are not k- 
corrected, which could well explain the relative redening of 
the faint objects. 



5.3 The local CfA2 Luminosity Function (LF) 

Zwicky magnitudes have been used to estimate the LF in 
the CfA2 redshift catalogue (Marzke et al. 1994). Marzke 
et al. used a Monte Carlo method to estimate how the LF 
parameters would change if the dispersion aM were 0.65 mag 
(closer to what we find here than the nominal 0.3 mag they 
used). They found that the true M« should be about 0.6 mag 
brighter, and that a similar conclusion would be reached 
for a small scale error. Marzke et al. used this estimate to 
conclude that a combination of incompleteness and a small 
(0.2 mag) scale error would be sufficient to move the CfA2- 
North values to those found from CfA2-South. 
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- m, < 15.5 




B-R 



Figure 16. Frequency distribution (in per cent) of galaxies as a 
function of CCD colour B — R for galaxies selected in different 
ways: (a) Zwicky sub-sample (top), (b) bright galaxies with CCD 
magnitudes B < 16 (middle) and (c) faint galaxies with CCD 
magnitude 19 < B < 20 (bottom). 



A detailed analysis of the implications of our new 
Zwicky magnitude calibration on the LF would require the 
redshift information and this is left for future work. It is nev- 
ertheless possible to use the mean relation that we found to 
show how we expect the LF to change. 

In practice, when fitting the Schechter parameters in 
equation 111 to observational data, the value of 0, is corre- 
lated with the values of M* and a. In our case we do not 
use a fit to data, but just model how the LF changes with 
an homogeneous linear shift in the magnitude scale. A zero 
point shift Zq in the magnitude scale, as in equation IItI will 
just change the value of Af* : 



Mt (corrected) = Af* + Zq 



(19) 



As the LF measures the number density of galaxies per mag- 
nitude inverval, a linear change in the magnitude scale will 
shift the amplitude of the LF proportionally to the shift in 
the magnitude interval, i.e. by A in equation (K&t). Thus we 
have: 



t (corrected) = A < 



(20) 



Thus the scale and zero point error in the Zwicky system 
give the following corrections for the CfA2 South results of 
Marzke et al. (1994) 

(fj^corrected) = A0. ~ 0.0124 ± 0.006 /i^ Mpc"^ (21) 

Af, (corrected) = M, + Zo ~ -19.3 ± 0.3. (22) 

If applied to the LF for the whole CfA2 survey: 

0. (corrected) = A </>. ~ 0.025 ± 0.009 h^ Mpc"^ (23) 

M^corrected) = Af. + ^o ^ -19.1 ± 0.2, (24) 

where we have added the errors in quadrature. Thus the 
corrected values are now closer the SAPM and LCRS (see 
Table pj). This can also be seem in Figure h/a, where we 
compare the corrected CfA2 luminosity function given above 
with the one corresponding to the SAPM. Note that the 





-22 



Figure 17. Luminosity function estimations. The continuous 
lines enclose the 2-sigma region in the SAPM estimation (bj 
band) whereas the dashed lines enclose the 2-sigma region in the 
CfA2 estimation (Bqcd) after correction for the model of the 
Zwicky magnitude system and errors. The upper panel shows the 
corrected LF for CfA2-South, and the lower for the whole CfA2 
survey. 



later is in the APM hj band, so that there could be some 
additional zero-point shifts between them (see §2). 



6 CONCLUSION 

In this paper we have presented CCD magnitudes for galax- 
ies around 204 Zwicky galaxies which sample over 400 square 
degrees, extending 6 hours in right ascension. This subsam- 
ple is drawn entirely from volume V of the Zwicky cat- 
alogue. We find evidence for a significant scale error as 
pointed out by Bothun & Schommer (1982) and Giovan- 
nelli & Haynes (1984). This is found by direct likelihood 
analysis, corrected for Malmquist bias, and also by notic- 
ing the angular correlations between the errors (Fig.y_0|) or 
the long tail of small objects in the frequency distributions 
of sizes (FigJld). The mean scale magnitude error is quite 
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large: Attiz — 0.62 ABccd, i-e. an error of 0.38 mag per 
magnitude, while the mean zero point is about —0.35 (equa- 
tion |l7| ), in agreement with Efstathiou et al. (1988). 

An illustration of the scale error in the Zwicky system 
can be seem in the galaxies shown in Figures La and h_l| 
These illustrations are strongly suggestive of an observer- 
bias effect, whereby some fainter galaxies are included in 
the sample due to their proximity to brighter objects, e.g. 
object #9 in Figure [13. This bias can also be seem in sta- 
tistical terms as a long tail of small objects in the frequency 
distribution of FigJlSl 

Huchra (1976) found only a 0.08 mag per magnitude 
scale error in a callibration of Zwicky galaxies with a pho- 
toelectric photometry of 181 sample of Markarian galax- 
ies, which are preferentially spirals. As spiral galaxies are 
bad tracers of clusters or groups, it is unlikely that these 
Markarian galaxies include any of the fainter galaxies that 
contribute to the proximity observer-bias mentioned above 
(which are mostly ellipticals). This effect would be hard to 
notice with photoelectric photometry which samples only 
one object at a time. Note that our analysis is restricted to 
a narrow range of magnitudes 13.5 < mz < 15.5, as com- 
pared to the wider range in Huchra (1976), but this narrow 
range contains the majority of the galaxies with mz < 15.5 
and therefore dominates all the relevant statistical proper- 
ties (such as the luminosity function). 

Bothun & Cornell (1990) have studied the calibration 
of the Zwicky magnitude scale using a sample of 107 spiral 
galaxies. They suggest that the errors in mz are minimized 
if mz is an isophotal magnitude a,t B — 26.0 mag arcsec"^, 
although it is clear from their Figure 2 that even within 
such a small sample of objects there is a 5 mag range of 
isophotes which give isophotal magnitudes corresponding to 
m.z- This suggests that our isophotal detection limit should 
be optimal for this comparison. 

Takamiya et al. (1995), also with photoelectric pho- 
tometry, find evidence for a large scale shift between Vol- 
ume I and Volumes II and V of the Zwicky catalogue. This 
shift appears to be of order 0.5 mag mag"'^ over the range 
14 < mz < 15.7, which is comparable to our findings for 
Volume V, although they find very little effect for Volumes 
II and V. However, we note that Figure 4 of Takamiya et 
al. shows B — mz vs. B, rather than mz- A similar rep- 
resentation of Figure of this paper looks very similar to 
Figure 4a of Takamiya et al., which is reasonable, given that 
Figure H implies a strong compression of the mz axis, which 
effectively hides the scale error. Unfortunately, the overlap 
between the 204 objects in our sample and the 155 objects 
in the Takamiya et al. data is too small to draw any de- 
tailed comparison, given that there are ~ 600 Zwicky galax- 
ies within the region of our survey. 

We have also estimated how this sc ale error could 
change the CfA2 luminosity function in §5.3. Figure |l7| 
shows how the corrected estimation is now closer to other 
local estimates, such as the SAPM. Finally, in 9a we give 
properties of the galaxies in our sample: colours, sizes and 
ellipticities, providing one of the largest samples of this kind. 
The local colour frequency distribution can be well approx- 
imated by a Gaussian distribution with mean B — R — 1.45 
and rms deviation of o" = 0.40. These colours compare well 
with synthetic B—R colours presented Fukugita et al (1995). 
But there is a significant number of galaxies with B — R> 2 



which are preferably found in clusters and groups, while 
most of the bluer objects with B — R < 1 (spirals to irreg- 
ulars) seem more isolated. These local properties are inter- 
esting in the context of galaxy evolution and star formation 
rates. This will be studied in more detail in future work. 
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